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highest.  The  simulations  show  that  such  effects  exhibit  a  pronounced  threshold  due  to  seif  focusing, 
and  disappear  when  the  focal  width  becomes  comparable  to  the  length  of  the  gain  medium.  In  the 
opposite  limit,  where  the  pump  radiation  at  the  lens  is  dominated  by  small  scale  transverse  structure, 
the  backscatter  can  reproduce  this  structure  in  detail  (i.e.,  exhibit  wavefront  reversal)  when  the  far 
field  of  the  pump  produces  an  interference  pattern  in  the  plasma.  The  plasma  then  behaves  as  an 
active  volume  hologram  as  it  amplifies  the  backward-propagating  random  noise  fields.  If  the  target  is 
moved  far  out  of  focus,  so  that  the  pump  components  can  no  longer  interfere,  the  ray  retracing  dis¬ 
appears,  and  the  actual  backscatter  profile  depends  upon  the  noise  structure.  Wavefront  reversal  is 
found  to  be  influenced  by  several  other  factors,  including  focal  spot  size,  spatial  gain  narrowing,  gain 
inhomogeneity,  and  self  focusing.  The  simulations  show  that  it  is  most  pronounced  when  the  spatial 
profiles  of  the  pump  and  the  backscatter  are  wide  in  comparison  to  the  transverse  wavelengths  of  the 
interference  structure,  and  the  1/e  gain  length  (especially  in  the  last  one  or  two  e-foldings)  is  large  in 
comparison  to  the  optical  diffraction  lengths.  It  is  enhanced  significantly  if  the  plasma  has  a  monoto- 
nically  decreasing  gain  coefficient  away  from  the  target.  In  the  strong  damping  limit,  where  the  non¬ 
linear  phase  change /cm  is  comparable  to  the  gain  coefficient,  filamentation  can  be  important  near  the 
target;  however,  its  effect  on  small  scale  ray  retracing  is  usually  minimal,  provided  that  the  above  con¬ 
ditions  are  satisfied.  If  the  focal  spot  decreases  below  a  critical  width,  however,  whole  beam  self 
focusing  takes  over  and  quickly  destroys  the  small  scale  retrace. 
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I.  INTRODUCTION 


Optical  ray  retracing,  in  which  the  transverse  spatial  pattern 
of  an  incident  beam  is  reproduced  by  its  backscatter,  is  a  well  known 


feature  of  stimulated  Brillouin  scattering.  The  effect  has  been 


|.ji  5~6 

studied  in  liquids  and  extended  plasmas,  as  well  as  in  laser-produced 


plasmas 


7-10 


In  the  latter  case,  it  occurs  even  when  the  critical  surface 


is  oriented  to  produce  a  strong  off-axis  specular  reflection,  as  shown 

10 


by  the  experiment  illustrated  in  Fig.  1.  This  paper  describes  numerical 
simulations  of  ray  retracing  in  laser-produced  plasmas,  using  a  steady 


11 


state  2-D  propagation  code  (BOUNCE).  In  its  present  form,  the  code 


12  13 

treats  the  strong  (ion-acoustic)  damping  limit,  ’  3  and  includes  self 


focusing  due  to  the  pondermotive  force; ^  however,  it  ignores  pump 

12 


depletion,  self  consistent  ion  heating,  and  ion  wave  saturation 
effects , 1 5”1 7 


It  is  useful  to  classify  ray  retrace  phenomena  roughly  into  two 
limiting  cases.  In  the  "whole  beam"  case,  the  pump  field  incident  at 
the  lens  has  a  broad  smoothly  varying  spatial  profile  (e.g.  Gaussian 
or  Hyper-Gauss i an)  that  can  be  focused  to  a  long  narrow  waist  within 
the  plasma.  The  narrow  waist  selectively  amplifies  only  those  components 
of  the  incident  noise  spectrum  that  are  propagating  back  along  the  near- 
axial  directions,  where  the  net  gain  is  highest.  The  simulations  show  that 
the  ray  retracing  associated  with  this  effect  can  exhibit  a  pronounced 
threshold  due  to  self  focusing,  and  disappears  when  the  focal  width 
becomes  comparable  to  the  length  of  the  gain  medium.  They  also  show 
that  although  the  angular  distribution  of  the  backscatter  at  the  lens 


is  qualitatively  similar  to  that  of  the  pump,  it  will  usually  be  somewhat 


18 


broader.  This  arises  from  spatial  gain  narrowing  of  the  backscatter 


Note:  Manuscript  lubmitted  January  24,  1980. 


profile  within  the  plasma,  which  leads  to  a  greater  diffraction  spreading 
as  the  radiation  returns  to  the  lens. 


In  the  opposite  limit,  where  the  pump  field  is  dominated  by  small 

scale  transverse  inhomogeneity,  the  backscatter  can  reproduce  this  structure 

in  detail  under  certain  conditions.  An  example  of  this  is  the  well- 

known  "question  mark"  experiment  performed  by  Eidmann  and  Sigel,^  in 

which  the  Brillouin  backscatter  from  a  laser  produced  plasma  retraces 

an  intensity  pattern  imposed  on  the  pump  beam  by  an  apertured  mask 

near  the  focusing  lens.  The  phenomenon  has  also  bein  observed  in  the 
k  1 9 

Brillouin  and  Raman  backscatter  from  liquids.  Ray  retracing  behavior 

in  which  the  spatial  structure  at  the  lens  resides  entirely  in  the 

1  2  20 

phase  rather  than  the  irradiance  has  been  studied  extensively,  ’  * 
and  the  feasibility  of  using  this  effect  to  correct  phase  aberrations 
in  a  laser  system  has  been  demonstrated.^ 

In  ideal  cases,  small  scale  ray  retracing  can  be  described  as 
wavefront  reversal  (or  phase  conjugation),  in  which  the  backscatter 

^  JL  i  ^ 

amplitude  (r)  is  proportional  to  the  conjugate  amplitude  BQ  (r)  of  the 

1  1  21-23 
incident  radiation.  Theoretical  studies  of  Zel'dovich  and  others  ’ 

indicate  that  phase  conjugation  occurs  under  conditions  where  any  small 

»'(  -f 

fraction  of  the  initial  noise  field  that  happens  to  be  proportional  to  £q  (r) 
can  grow  faster  than  the  remaining  portions,  and  thus  eventually  predominates. 
The  physical  basis  for  wavefront  reversal  can  be  readily  visualized  in 
the  case  where  the  pump  profile  at  the  lens  consists  of  two  or  more 
isolated  hot  spots.  The  beams  from  these  hot  spots  overlap  near  focus 
to  create  an  interference  pattern  in  the  plasma.  This  pattern  selectively 
amplifies  only  those  noise  components  propagating  oppositely  to  any  one 
of  the  pump  beams,  and  simultaneously  Bragg-d; ff racts  energy  back  along 


11  24-26 

all  of  the  remaining  beams.  '  If  the  plasma  is  moved  far  out  of 

focus,  so  that  the  beams  no  longer  overlap,  then  the  wavefront  reversal 
disappears,  and  the  actual  backscatter  profile  depends  upon  the  incident 
noise  structure. 

To  examine  this  mechanism  in  greater  detail,  consider  the  simple  model 
shown  in  Fig.  2.  A  pair  of  isolated  pump  beams  with  propagation  vectors 
k  and  k  ,  and  initial  widths  A~  are  focused  into  a  plasma  slab 

O  O  D 

of  thickness  L>:>A  .  In  the  focal  region,  the  beam  waists  (indicated 

y  2 

by  the  dashed  lines)  overlap  to  produce  an  interference  pattern  |£>(r)| 

of  spatial  frequency  K=  kQ  -  kQ,  width  d  —  2*of/Ag>>AQ  an^  length 

8\f  /A  >L,  where  f  is  the  focal  length  of  the  lens.  A  plane  wave 
O  D 

-f 

noise  field  incident  along  some  backward  direction  k^  therefore  sees 

a  periodic  gain  coefficient  proportional  to  |  £q(?)  1^  >  as  indicated 

by  the  parallel  lines  in  Fig.  2.  If  the  dimensions  L  and  d  are  large 

in  comparison  to  the  periodicity  A*  2VK,  then  the  region  behaves  as 

27 

an  active  volume  hologram.  The  k^  wave  will  grow  most  rapidly,  while 

-f  I 

simultaneously  Bragg-d i f f ract i ng  energy  into  the  mode,  when  this 
pair  of  coupled  waves  satisfies  the  Bragg  condition 

k,  -  f,'  k0-  K  (.) 

28 

within  the  medium.  The  gain  is  maximum  when  the  periodic  intensity 


I  c  /  1 2  -*■  -►  *  t 

pattern  |Cj(r)|  developed  by  the  k^ ,  k^  waves  is  exactly  in  phase 

with  that  of  I  t  (?)  I ^ :  i . e. ,  when  k.  =  -k  and  k  '  =  -k  1 .  Si  nee 

kj 1  =  kj-  kQ  for  Brillouin  scattering,  this  condition  occurs  when 

0^  and  ©1  'ac90'  -  Any  arbitrary  phase  difference  <t>  between  the 

pump  beams  will  also  be  reproduced  in  the  two  backscattered  beams 

(i.e.,  £j(r)~£o’'  (r))  because  the  backscatter  phase  difference  depends 

upon  the  phase  of  the  interference  pattern,  which  in  turn  depends  upon  <(>. 


3 


A  coupled  wave  analysis 


has  shown  that  the  intensities  of 


11 ,21* 

the  two  backscattered  beams  become  nearly  equal  after  a  few  e-foldings, 
and  that  the  gain  maximizes  at  8^  ■  8g  (assuming  ■  kg)  with  a  1/e  angular 

bandpass  68 ^  ~  's  defined  in  Fig.  2).  For  a  uniform 

medium  of  thickness  L  and  average  amplitude  gain  coefficient  a>>L  \ 
one  obtains  68j  ■  (8/aL)^  a/kQ  ®g«  giving  an  angular  resolution 

(8q-+68i ) 

<<h)>.|„  -  W«U*  (2.) 

For  a  semi-infinite  medium  of  exponentially  decreasing  gain  a(z)  = 
aQ  exp  (-z/ji),  where  a(J£  >>  1,  one  obtains 

S>.l„  ■  <2'»'ko*>*  •  <2b> 

As  long  as  the  widths  of  the  focal  spot  and  backscatter  remain  large 

in  comparison  to  the  transverse  period  A,  the  wavefront  reversal  will 

be  evident  if  ®0  >:>(0D)mjn.  F°r  a  fixed  gain,  Eq.  (2a)  requires  that 

2 

the  diffraction  length  1/kQ  8^  must  be  much  shorter  than  the  effective 

21  22  ■+• 
gain  length  1/a,  ’  or  in  terms  of  the  spatial  frequency  K  of  the 

interference  pattern, 

K2/kQ  »  a.  (3) 

The  physical  argument  for  this  requirement  can  be  stated  as  follows:  As  long 
as  there  is  enough  gain,  the  transverse  intensity  variations  of  the 

pump  will  always  tend  to  "pull"  the  backscatter  intensity  into  a  similar 

.  22 
pattern  as  it  propagates  back  through  the  active  medium.  However,  if  this 

pulling  effect  occurs  in  too  short  a  distance,  the  backscatter  phase 

(which  cannot  be  directly  affected)  will  remain  similar  to  that  of  the 

incident  noise,  which  is  uncorrelated  with  the  pump  phase.  The  backscatter 

would  then  begin  to  diverge  from  the  pump  as  soon  as  it  leaves  the  active 

medium.  (Within  the  context  of  the  above  model,  the  medium  would  behave 

4 


mm 


as  a  simple  diffraction  grating  with  many  orders,  rather  than  a  volume 
hologram.)  Now  consider  the  case  where  the  gain  occurs  slowly  in 
comparison  to  diffraction.  The  amplification  then  takes  place  over 
a  long  path  that  effectively  filters  out  all  of  the  backward 
propagating  modes  excepting  the  one  that  would  have  developed  the  transverse 
intensity  pattern  even  without  the  pulling  effect;  i.e.,  the  mode 
proportional  to  (r). 

The  numerical  calculations  presented  here  will  illustrate  the 

main  results  of  the  volume  hologram  model,  and  explore  the  influence 

.  11  18,29 

of  several  factors  (e.g.  focal  spot  size,  spatial  gain  narrowing, 

longitudinal  gain  inhomogeneity,  and  fi lamentation)  on  the  small  scale 

ray  retracing.  The  simulations  show  that  small  scale  retracing  is 

most  pronounced  when  (i)  the  power  gain  is  >1000,  (ii)  the  focal  spot 

and  backscatter  remain  wide  in  comparison  to  the  period  of  the  transverse 

interference  structure,  and  (iii)  the  1/e  gain  length  (especially 

in  the  last  one  or  two  e-foldings)  is  large  in  comparison  to  the  diffraction 

length  of  the  interference  structure.  It  is  enhanced,  for  example, 

if  the  plasma  has  a  monotonical ly  decreasing  density  (and  therefore  gain 

1 1 ,2k 

coefficient)  away  from  the  target.  In  the  strong  damping  limit,  where 

the  nonlinear  phase  change/cm  is  comparable  to  the  gain  coefficient,^ 

14  31 

fi lamentation  can  be  important  near  the  target;  '  however,  its 
effect  on  small  scale  ray  retracing  is  usually  minimal  (with  the  density 
gradient)  if  the  above  conditions  are  satisfied.  If  the  focal  spot 
decreases  below  a  critical  width,  however,  whole  beam  self  focusing  takes 
over  and  quickly  destroys  the  small  scale  retrace. 

In  the  following  section,  the  reduced  wave  equations  are  derived 
for  the  steady  state  pump  and  backscatter  amplitudes  in  the  strongly 
damped  limit,  assuming  paraxial  ray  conditions.  Section  III  describes 


5 


the  propagation  code,  and  compares  it  to  an  analytic  calculation  for 
a  simple  test  problem.  Section  IV  presents  numerical  simulations 
with  1.06m  radiation  for  several  cases,  including  both  whole  beam  and 
small  scale  retrace  effects.  These  results  are  summarized  In  Sec.  V. 
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I I .  THEORY 


Stimulated  Brillouin  scattering  arises  from  the  interaction  between 

a  pair  of  optical  waves  Eg(r,t),  Ej(r,t)  of  frequencies  and  an 

ion-acoustic  density  fluctuation  <$n(r,t)  of  frequency  -  u>j .  This 

interaction  can  be  described  by  Maxwell's  equations  and  the  two-fluid 

plasma  equations  (neglecting  electron  inertia)  with  an  oscillating  ponder- 
17  32  33 

motive  force.  Similarly,  intensity-dependent  variations  of  the 

time-averaged  plasma  density  n(r)  can  be  related  to  the  quasi-static 
pondermotive  pressure.  For  |<5n(r,  t)|«  n(r),  the  results  are 


32 

(  —  + 
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+  Aw 
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cs2  72)  <Sn 

[  6n  V(E2  +  E2)av  +  2nV(EQ  E,)], 


(4a) 


(4b) 


where  Aws  is  the  acoustic  energy  damping  rate, 
cs  =  [(ZT,  +  3T.)/M]* 


(5) 


is  the  acoustic  propagation  velocity,  Z  is  the  ion  charge,  m(M)  is 

the  electron  (ion)  mass,  Tg(T.)  is  the  electron  (ion)  temperature, 

and  (  )  denotes  a  time  average  over  one  optical  cycle.  The 

first  two  terms  on  the  right  hand  side,  which  arise  from  the  quasi-static 

pondermotive  forces,  are  usually  ignored  in  1-D  treatments  where 

the  spatial  variation  of  the  pump  intensity  is  sufficiently  weak. 

If  the  acoustic  transit  time  across  the  focal  spot  is  short  in  com¬ 
parison  to  the  pulsewidth,  then  n(?)  is  determined  by  the  quasi-static 


7 


pressure  balance: 


14 


V{(Te  +  T./Z)n} 


+  (n/8™c)  V  (^2  +  E2)av  =  0, 


(6) 


2  2 

where  nc  =  mu)o  /4ire  is  the  critical  density, 
the  solution  is 

r  <Eo  +  E?>*v  i 
n  r)  =  n  (r)  exp  [  -  g  {j  -f.'/z)'  J 

c  e  i 


For  isothermal  conditions, 


(7) 


where  n  (?)  is  the  static  density  at  low  irradiance. 
o 

We  now  assume  paraxial  ray  conditions;  i.e.,  allow  the  pump  (Eg) 
and  acoustic  (Sn)  waves  to  propagate  approximately  along  the  -z 
direction,  while  the  backscatter  (E^)  travels  approximately  along  +z. 

For  simplicity,  the  condition  nQ  =  nQ(z)  w'  1 1  also  be  assumed.  Following 

34 

Kauffman  and  Cohen,  one  can  then  remove  the  rapidly  oscillating  portions 
of  Eqs  (4)  by  expanding  in  terms  of  the  complex  wave  amplitudes  £Q,  £1 
and  6n: 


E.(r,t)  >  i  6.  (?,t)  exp  {  i o.  J  k. 
J  J  J  J 


(z 1 )  dz 1  -  i  uk  t  }  +  c.c.  (8a) 


6n(r,t)  =  i  6n(r,t)  exp  { 


z 

-i  J  k(z')  dz1  -  i  to  t  }  +  c.c.,  (8b) 


where  aQ=  -1,  =  +1 ,  and  the  carrier  frequencies  and  propagation 

vectors  satisfy  the  relations 


u  =  <»>0  -  o>j  ,  k(z)  =  kQ(z)  +  k1  U)  —  2kg (z) 


c2  k?  (z)  =  u)2  -  4ire2  n  (z)/m. 
J  J  o 


(9a, b) 
(10) 
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■ 


The  acoustic  amplitude  (and  thus  the  backscatter)  cannot  grow  significantly 
unless  k(z)  satisfies  the  resonance  condition^ 

k(z)  ks  =  ui/cs;  (11) 

36 

i.e.,  the  k-vector  mismatch 

Ak(z)  =  k (z)  -  ks  (12) 

should  satisfy  |a  k(z)  |«  ks 

Expressions  (8) - ( 1 2)  are  now  substituted  into  Eqs  (4a,  b) , 
subject  to  the  restrictions 

kQ  >>a*  K2/k0  >>  -  (13) 

where  is  the  gradient  scale  lengths  of  k(z).  The  equations  for 
the  (slowly-varying)  steady  state  amplitudes  then  reduce  to 


i  „  2 , 
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where  ^  =(3/3x,  3/ 3y)  is  the  transverse  gradient  operator,  V j  ( r )  h 
e  £j(r)/mw.  are  the  electron  quiver  velocity  amplitudes, 


9 


v  =  {T  /m)5  Is  the  electron  thermal  velocity,  e  h  1  -  n  /n  the 
e  e  o  c 

plasma  dielectric  constant,  and  n=  1  +  3  Tj/ZT  .  (On  the  right  hand 

sides  of  Eqs  (Ik),  we  have  approximated  kj  —  kQ=  o)q/c  and  k  c£  k^  . ) 

2 

The  V.J.  terms  in  these  equations  account  for  diffraction  effects,  while 
the  (nQ  -  n)  terms  describe  the  intensity-dependent  phase  shifts  and  self 
focusing.  This  intensity  dependence  is  given  by  Eq.  (7),  which  can 


be  written  as 


-  -  r  1 V  |2  +  |v  I2, 

n (r)  =  Mr)  exp - = -  1  , 

°  L  kn1  v  2  J 


where  n  =1  +  T . /ZTg. 

To  evaluate  the  final  terms  of  (lka,b),  one  must  solve  (Ike)  for 

12  13 

6n(r) .  The  present  paper  will  deal  only  with  the  strong  damping  limit  ’ 
where  Atos  predominates  over  the  convective,  diffraction,  and  phase 
mismatch  terms  in  (ike);  i.e., 


Aw  /c  »  a,  K  /k  ,  Ak. 
s  s  o 


6n(r)/n(r)  Si 


-i  (Q/2H)  Vq  V^'/v^, 


where  Q  =  •%  cs^Aws  =  “/Au)s  can  be  interpreted  as  the  "quality  factor"  of 
the  ion-acoustic  resonance.  Condition  (16)  is  most  likely  to  be 
found  in  the  low  density  corona  region  where  ion  Landau  damping  is 
heavy,12  and  the  plasma  scale  lengths  are  large. 1(^  The  effects  of  the  phase 


mismatch  on  optical  ray  retracing  will  be  examined  in  a  later  publication. 


Substituting  expression  (17)  into  Eqs  (14a,  b) ,  one  finally  obtains 


7?>  Vo 


n  -n 
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V  + 
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(18a) 
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v2)  V 
V  V1 


n  -n  k  Q 

-2—  v  +  _2_ 
n£  1  8en 


(18b) 


e 

The  final  terms  of  (l8a)  and  (18b)  describe  the  pump  depletion  and 

gain,  respectively.  Excepting  for  these  (and  the  small  difference  between 

kj  and  kQ)  the  equations  for  and  Vq  would  be  identical.  This  is  the 

reason  why  the  backscatter  field  amplitude  £1  =  (m/e)^  can  become 

proportional  to  the  conjugate  pump  amplitude  £' q  =  (m/e)u>  when 

the  gain  term  is  small  in  comparison  to  the  diffraction.' 

In  solving  Eqs.  (l8),  some  care  must  be  taken  to  ensure  that 

the  condition  | 6 n | < <n  remains  satisfied;  i.e.,  according  to  expression 

(17),  the  formalism  cannot  simultaneously  treat  heavy  pump  depletion 

2  2 

(  |v ^  |»|vj)  and  large  static  pondermot i ve  density  variations  (  |v|  «  J*vg) . 
The  calculations  presented  in  this  paper  will  ignore  pump  depletion 
entirely,  and  thus  assume  | V ^ | << | VQ I  .  Since  most  of  the  examples  chosen 
here  will  also  satisfy  the  condition  I Vq| 2<<  4v2,  it  is  instructive  to 
write  Eqs  (18)  with  these  approximations: 
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V1 

n 

n 

c 

2  V1 

V 

e 

(19b),  one 

then 

sees 

that 

the 

ratio 

of  the 

gain  to  self  focusing 

terms  becomes  (n'/n)  Q~  Q  >1 ,  as  noted  in  Sec.  I  and  Ref.  30. 
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A  complete  solution  of  Eqs  (18)  or  (19)  would  require  an  absorption 
and  hydrodynamic  model  that  relates  noCr),  T  (r)  and  Q(r)  to  the 

12 

incident  irradiance,  and  would  Include  self  consistent  ion  heating 
and  various  ion  wave  saturation  mechan? sms . **’ ' ^  ^  Initially,  however, 
we  wish  to  concentrate  on  the  role  played  by  diffraction  and  self  focusing 
in  the  ray  retracing.  The  examples  presented  here  will  therefore  deal 
with  the  simplest  model  where  nQ,  T&  and  Q  are  either  constants  or  pre¬ 
determined  functions  of  position.  The  inclusion  of  a  self  consistent 
hydrodynamic  model  or  ion  wave  saturation  must  be  reserved  for  a  later 


III.  PROPAGATION  CODE 


In  its  present  form,  BOUNCE  solves  Eqs  (18)  (ignoring  the 
2  2 

I  V,  I  /v  term  in  18a)  in  two  cartesian  coordinates  (x,z),  using  a 
1  l  1  e 

fast  Fourier  transform  subroutine^  to  handle  the  transverse  (x)  variation. 
If  Vj (* ,z)  is  the  Fourier  transform  of  Vj (x,z) ;  i .e. , 


Vj(k,z)  =  T  {Vj (x,z) }  ,  Vj(x,z)  =  T  1  (V.(k,z)}  , 


(20) 


then  Eqs  ( 1 8 )  become 
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where  n  is  given  by  Eq.  ( 1 5) - 

As  we  have  shown,  the  wavefront  reversal  is  most  effective  when 
2 

the  diffraction  term  k  /2k^  is  large  in  comparison  to  the  gain  and 
self  focusing  terms.  This  rapidly  varying  part  can  be  effectively 
removed  from  the  integration  by  using  the  unitary  transformations 


Vj(k,z)  =  V j 1 (< , z)  exp  (-i0j  x  z/2kj), 


(22) 


(where  =  -1  and  a,  =  +1)  to  rewrite  (21a, b)  in  the  "interaction 
representation" 


9V 
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9V, '  .2  ^  ik  n  -n  k  Q 
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—  =  exp  (~2k7  }  T  {  TF  “n  V1  ^  ~ 

1  c  c 


o I  V,  }  >  (23b) 
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These  equations  are  then  integrated  by  a  predictor-corrector  technique, 
using  Eqs  (20)  and  (22)  to  obtain  Vo(x,z)  and  Vj(x,z)  at  each  new 
point  aiong  z. 

The  input  parameters  include  the  density  nQ(z),  quality  factor  Q, 

temperature  ratio  T./T  wavelengths  and  incident  optical  field 

1  e  j 

amplitudes  £.(r)  -  (mo^/e)  Vj(r).  (See  Fig.  3)  In  all  cases  treated 
here,  we  have  chosen  Xq  -  X^  ■  1.064y,  and  for  simplicity,  have 
appoximated  e  «  1  and  n  =  n'“  1.  The  pump  field,  whose  transverse 
profile  is  usually  defined  at  the  lens  plane,  is  normalized  to  a 
specified  average  ratio  (|VQ|  )  ^/v  *  at  t*ie  backscatter  exit  plane 

z2  of  the  active  medium.  The  average  of  any  quantity  ^(x.z)  related 
to  the  amplitude  Vj(x,z)  is  defined  as 

<^>2  5/*j(x,z)  |V.(x,z)[2  d)//|Vj(x,z)|  2dx,  (24) 

and  the  axial  position  z  is  measured  from  the  focal  plane  of  the  lens. 

The  backscatter  "noise"  field  C^x.zl),  which  is  arbitrarily 

normalized,  is  specified  at  the  entrance  plane  zl .  Its  x  dependence 

was  in  most  cases  chosen  to  be  a  stochastic  function  with  a  random 

38 

phase  Gaussian  power  spectrum  of  (l/e)  angular  width  A6^.  For  most  of 
the  examples  discussed  here,  A0^  =  .25  radian;  however,  the  backscatter 
profiles  do  not  depend  strongly  on  A0^  for  power  gains  >100. 

The  code  plots  the  pump  and  backscatter  profiles  at  the  zl ,  z2  and 

lens  planes,  and  calculates  several  additional  parameters.  These  include 

2  i 

the  RMS  widths  2  (xj  )  *  of  the  pump  and  backscatter  beams,  an 
effective  power  gain 
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and  an  average  nonlinear  phase  shift 


In  most  cases,  G  and  B  satisfy  the  approximate 

as  one  might  expect  from  Eqs.  (19).  Finally, 

22  29 

pump-backscatter  correlation  function  * 


n  <|V  |2>  .  , 

o  ' '  o1  'z  dz.  ( 

£  n  2 

c  v 

e 

relationship  Jin  GCtf2QB, 
the  code  calculates  the 


|/£0(x)  £1(x)dx|2 

/le0(x)|2d>«J'  |S,W| 


where  the  field  amplitudes  are  evaluated  at  the  lens.  This  is  a  measure 
of  the  phase  as  well  as  amplitude  correlation  between  the  beams,  and 
in  the  limit  of  exact  wavefront  reversal  £^(x)~6^  (x) ,  one  obtains 
H  =  1. 

In  order  to  evaluate  the  accuracy  of  the  code  (and  provide  a 

graphic  illustration  of  the  volume  hologram  model)  we  compared  it 

to  the  analytic  theory  in  a  simple  test  problem.  A  uniform  slab  of 

thickness  L  =  100W  and  plasma  density  n  =  . 1  n  is  illuminated 

o  c 

directly  at  z2  =  +50U  by  a  pair  of  intersecting  plane  waves  with 
hyper  Gaussian  profiles;  i.e., 


€-q(x,z2)  =  A(0,z)  +  A  (0p, x)  exp  (  -ikQ  8^  x) 


(28) 


where  0D  ■  .16  radian, 


A(0,x)~exp  {  -  [2(x  -  i  L9)/d]  20  }  , 


(29) 


w 


and  the  width  d  is  l80u.  The  resulting  interference  pattern  (Fig.  4) 

1 1  29 

is  sufficiently  broad  and  flat  that  edge  effects  and  gair  narrowing  * 

are  unimportant  in  this  case.  The  quality  factor  Q  »1  and  average 

intensity  '(|V  I  >  „<<  v  c  were  chosen  to  given  an  effective  (amplitude) 

1  o 1  ' zZ  e 

gain  coefficient  a  ■  .0376m  1  with  negligible  self  focusing  (B  ■  0). 

The  backward  noise  field  for  this  case  was  a  hyper  Gaussian  plane  wave 

£j(x,z1  *  -50y)  =  A(01,x)  exp  (+ikQ  0^x)  (30) 

of  width  240y. 

For  01  =  0  (Fig.  4)  the  backscatter  at  z2  has  developed  a  periodicity 

exactly  in  phase  with  that  of  the  pump,  as  described  in  Sec.  I.  This 

periodicity  produces  the  peaks  around  6=0  and  .16  radian  back  at  the 

lens  plane,  as  one  would  expect  from  Eq.  (28).  For  0^  >  0,  the  backscatter 

pattern  at  z2  shifts  slightly  to  the  right,  thereby  reducing  the  overlap 

with  the  pump  and  lowering  the  gain.  '  At  the  lens  plane,  both  peaks  shift 

a  distance  01  to  the  right.  An  analogous  shift  and  lowering  of  gain  occurs 

when  0,  <  0.  The  resulting  tuning  curve  (Fig.  5)  shows  good  agreement  with 

1 1  24 

the  analytic  theory,  ’  and  illustrates  the  limitations  of  the  volume 
hologram  mechanism.  The  small  discrepancy  between  the  .08  radian  (1/e) 
width  of  this  curve  and  the  value  50tf(8/aL)^  a/kQ0D  =  .058  radian 
from  Eq.  (2a)  arises  from  the  fact  that  aL  =  3-67,  while  (2a)  is  strictly 
applicable  only  if  aL  »  1. 
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IV.  RESULTS 


In  all  of  the  remaining  examples,  the  pump  field  incident  at  the 

lens  is  assumed  to  be  uniphase.  Simulations  of  wavefront  conjugation 

experiments  where  the  spatial  structure  of  the  pump  resides  in  the 
1-3  20 

phase  *  will  be  addressed  in  a  later  publication. 

A.  Whole  Beam  Effects 

Whole  beam  ray  retracing  is  illustrated  in  Figs.  6-8.  In  Fig.  6a, 

the  incident  pump  profile  is  a  half  blocked  hyper-Gauss ian  function 

of  e  angular  width  .19  radian.  This  is  focused  into  a  300p  thick  plasma 

of  exponentially  decreasing  density  nQ  =  .2  nc  exp  (-z/£),  where  £  »  100P, 

with  the  backscatter  noise  field  incident  at  the  focal  plane  (zl  »  0) . 

To  ensure  negligible  self  focusing  in  this  example,  we  chose  (|  VQ|  )« 

2 

«ve  ,  but  allowed  Q  »1  to  obtain  a  high  power  gain  In  G  »  11.  The  long 

narrow  focal  waist  within  the  plasma  (d2?12P  at  zl  and  70u  at  z2  *  300u) 

amplifies  only  those  noise  components  travelling  back  along  the  pump  beam, 

as  seen  by  comparing  the  profiles  at  z2  and  the  lens.  The  backscatter 

profile  at  z2  also  shows  gain  narrowing  due  to  higher  amplification 

]  g 

in  the  more  intense  central  portion  of  the  pump  beam.  This  also 

occurs  in  all  subsequent  examples,  and  produces  the  Gaussian-1 ike  tails  back  at 

the  lens.  In  spite  of  such  effects,  one  obtains  good  field  correlation 

(H  =  .72),  with  only  5%  of  the  backscatter  hitting  the  blocked  portion 

of  the  lens,  in  qualitative  agreement  with  the  experiment  of  Fig.  1.  ® 

Similar  results  are  found  at  higher  gain,  with  somewhat  greater  broadening 
at  the  lens  due  to  more  gain  narrowing  at  z2. 
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Conditions  in  Fig.  6b  are  similar  to  those  in  the  above  example,  ex- 

2  2 

cepting  that  we  now  choose  Q  ■  1  and  (  jVQ|  )Z2^e  "  *^4  to  allow  significant 

self  focusing  (B  =  7.6)  along  with  the  high  gain  Jin  G  -  3.9.  This  is 
evident  in  the  fi lamentation  of  the  pump  profile  at  zl  and  some  additional 
narrowing  of  the  backscatter  at  z2;  however,  the  backscatter  at  the 
lens  remains  well  correlated  with  the  pump  (H  =  .70),  with  still  only 
5%  in  the  shadow  region. 

When  the  pump  field  is  focused  at  z2  =  0,  one  obtains  similar  ray 
retracing  behavior,  but  the  self  focusing  can  cause  a  pronounced  threshold 
effect.  An  example  of  this  is  shown  in  Figs.  7a,  b.  A  hyper-Gauss ian 

-  If 

pump  field  of  e  angular  width  .12  radian  is  focused  to  a  20p  central 

2  2 

spot  of  average  i ntensi ty ()Vo |  )z2^ve  =  -508  at  z2,  then  propagates 

through  a  400y  plasma  slab  of  uniform  density  nQ  and  Q  =  1.  For 

nQ  =  .01  n c ( F i g .  7a),  the  pump  beam  expands  to  about  45 y  at  zl,  with  self 

focusing  effects  playing  only  a  minor  role  (B  *  1.2).  Although  the 

geometrical  conditions  are  ideal  for  whole  beam  ray  retracing,  the 

power  gain  &n  G  =  1.3  is  too  low  to  provide  sufficient  spatial  mode  filtering 

of  the  broadband  noise  source;  consequently,  the  effect  is  minimal. 

In  Fig.  7b,  the  density  has  been  doubled  to  nQ  =  .02  nc>  increasing 
the  self  focusing  (B  =  2.9)  and  reducing  the  central  beam  diameter  at 
zl  to  about  30u.  The  resulting  increase  in  pump  intensity  within  the 
slab  raises  the  gain  to  in  G  =  7.3,  and  produces  the  clear  retrace  at 
z2  and  the  lens.  The  probable  reason  why  £n  G  increases  faster  than  B 
is  that  after  the  first  one  or  two  e-foldings,  the  backscatter  is  confined 
primarily  to  near-axial  backward  modes,  which  are  then  amplified  more 
effectively  in  the  remaining  active  material.) 
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Figure  8  Illustrates  a  case  where  ray  retrace  would  not  be 
expected  to  occur.  Here,  the  plasma  is  a  uniform  slab  of  density  n  •  .1  nc, 
which  lies  between  300  and  **00 y.  The  thickness  is  smaller  than  the 
l80-2**0u  (e  )  beam  width  of  the  pump,  and  cannot  provide  effective  filtering 
action  to  favor  the  backward  direction.  Thus,  in  spite  of  high  gain  tn  G  «  7- 
the  ray  retrace  is  minimal,  and  the  backscatter  seen  at  z2  is  essentially 
an  amplified  version  of  the  noise  field  (whose  angular  width  was  broadened 
to  .5  radian  for  this  example).  Additional  evidence  for  the  amplified 
noise  character  of  the  scattered  radiation  is  the  low  correlation  Ha.OOl. 

B.  Small  Scale  Effects 

In  the  following  examples,  the  pump  angular  intensity  dis- 
tribution  at  the  lens  consists  of  three  isolated  Gaussian  peaks  of  e 
width  Ag/f  =  .025  radian,  with  the  maxima  at  -.066*1,  and  .0117  and  .121 
radian.  The  nearest  neighbor  spacing  is  thus  .078  radian,  while  the  total 
angular  spread  is  .187  radian. 

In  Fig.  9,  the  pump  radiation  is  focused  into  a  uniform  lOOu  slab 

of  plasma  density  nQ  =  . 1  nc>  with  a  quality  factor  Q>>1  and  average 
2  2 

intensity^) VQ|  )z2<<ve  chosen  to  produce  high  gain  *,n  G  ■  12.5  but 
negligible  self  focusing  (B  -  0) .  The  three  interfering  pump  beams 
at  the  focal  plane  zl  ■  0  produce  a  quas i-per iodic  intensity  pattern, 
which  is  superimposed  upon  a  Gaussian  envelope  of  total  width  d<s2Xof/Agae85p . 

At  the  z2  plane,  this  envelope  has  expanded  and  shifted  somewhat  due 
to  the  small  divergence  and  separation  of  the  focal  waists.  Near  the 
axis,  where  the  average  gain  is  high,  the  backscatter  has  begun  to  develop 
an  intensity  pattern  in  phase  with  that  of  the  pump,  just  as  it  did 
in  the  test  problem  (Fig.  **) .  In  this  case,  however,  the  spatial  gain 
narrowing  ’  of  the  backscatter  envelope  is  more  severe,  giving  a 
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total  envelope  width  of  only  as  30p. 

The  ray  retracing  is  readily  apparent  back  at  the  lens  plane,  but 

the  maxima  are  slightly  displaced,  and  three  additional  peaks  of 

significant  height  are  also  present.  This  type  of  degradation  occurs 

because  the  angular  bandpass  of  the  gain  medium  (due  to  its  finite  thickness) 

is  too  broad  to  effectively  filter  out  all  of  the  non-retracing  components 

in  the  incident  noise  spectrum.  Using  Eq.  (2a)  with  the  average  amplitude 

2  2  “1 

gain  coefficient  a*(k  /8)  (n  /n  )  Q(|V  |  )  ,/v  =  .052**  p  ,  one  obtains 

O  O  C  O  ZZ  6 

(eJ  .  ®  . 1  radian,  which  is  comparable  to  the  angular  distance  between 

adjacent  peaks  in  the  pump  beam.  The  broadening  of  the  backscatter  peaks 

arises  from  the  gain  narrowing^  noted  in  the  preceeding  paragraph;  i.e., 

for  an  envelope  of  width  dj2£  30p  at  z2,  one  obtains  a  total  angular 

width  2AQ/d1c;  .07  radian  back  at  the  lens.  The  gain  narrowing  filters 

out  most  of  the  higher  spatial  frequency  structure (from  the  incident 

noise)  that  would  otherwise  get  to  the  lens  plane  within  the  finite 

angular  bandpass  of  the  gain  medium. 

Figure  10a  illustrates  an  enhancement  of  ray  retracing  when  the  uniform 

slab  geometry  is  replaced  by  a  density  gradient  of  equal  scale  length 

and  total  thickness  300y;  i.e.,  nQ  =  .2  n£  exp  (-z/J l) ,  where  i.  =  100y, 

2  2 

and  0^  z  <  300p.  We  again  choose  Q>>  1  and(|VQ|  )z2/,ye  <<:  ^  to  obtain 
in  G  =  9.2  with  negligible  self  focusing  (B  =  0).  For  this  case,  one 
can  apply  Eq.  (2b)  to  estimate  ( ©q) m ; - 063  radian,  which  should  allow 
the  three  peaks  to  be  fully  resolved  in  the  backscatter.  This  expectation 
is  confirmed  by  the  profiles  in  the  lens  plane.  The  three  amplitudes 
are  now  nearly  equal,  and  the  mismatch  and  additional  structure  evident 
in  Fig.  9  have  almost  completely  disappeared.  An  improvement  in  ray 
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retracing  is  also  evident  in  the  better  pump-backscatter  synchron'zatlon 
at  z2,  and  the  larger  correlation  number  H  *  .50,  compared  to  H  *  .33 
for  the  preceeding  case. 

The  physical  reason  for  such  enhancement  is  that  although  the 

plasma  beyond  100u  accounts  for  only  a  few  e-foldings  of  amplification, 

the  local  gain  coefficient  satisfies  condition  (3).  so  the  region 

can  still  provide  significant  transverse  mode  discrimination.  For  example, 

the  highest  (amplitude)  gain  coefficient  in  the  last  two  e-foldings  of 

intensity  (the  region  between  135  and  300v» )  is  .013w  ' ,  giving  a  1/e 

gain  length  of  77u.  From  the  known  angular  separations  at  the  lens, 

one  obtains  K  .  /k  *  .078  radian,  giving  a  maximum  diffraction  length 
min  o 

of  k  /K *  -  28u. 

o  min 

A  comparison  between  the  z2  profiles  in  Figs.  9  and  10a  reveals 

that  gain  narrowing  is  less  evident  in  the  gradient  geometry,  presumably 

because  of  a  greater  opportunity  for  diffraction  spreading  in  the  100-300P 

region.  This  results  in  slightly  narrower  backscatter  peaks  at  the  lens; 

however,  the  residual  broadening  is  still  evident,  and  it  appears  to  be 

the  primary  reason  why  H  <1.  similar  reduction  of  wavefront  reversal 

due  to  spatial  gain  narrowing  was  predicted  in  Ref.  29- 

In  Fig.  10b,  the  conditions  are  similar  to  those  of  Fig.  10a,  excepting 

that  we  now  choose  Q  =  1  and{|v  |^)  „/v  ^  *  .213  to  allow  self  focusing 

o  zt  e 

(B  =  5*8)  along  with  the  high  gain  tn  G  ■  10.0.  This  causes  some  fi lamentation 
and  shifting  of  the  pump  intensity  spikes  at  zl,  but  little  narrowing 
of  the  envelope.  The  effect  on  the  backscatter  profile  at  z2  is 
minimal.  Excepting  for  small  differences  in  the  amplitude  of  the  spikes 
and  a  slight  envelope  narrowing,  this  intensity  distribution  is  essentially 
the  same  as  that  of  Fig.  10a.  Back  at  the  lens,  the  ray  retrace  appears 
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to  have  actually  been  enhanced  somewhat,  and  the  correlation  has  increased 
to  H  =  .55. 


The  relatively  minor  influence  of  self  focusing  in  this  example 
can  be  qualitatively  understood  in  terms  of  Bespalov-Talanov  perturbation 
theory.  According  to  this  theory,  a  small  intensity  ripple  of  spatial 
frequency  k  grows  exponentially  at  the  rate 


/I  dB 

aSF  (k  dz 
o 


(3D 


o 

due  to  self  focusing  effects.  If  <  exceeds  the  cutoff  frequency 
k  =  (A  kQ  dB/dz)^,  then  becomes  imaginary  and  the  growth  ceases; 
i.e.,  the  diffraction  spreading  begins  to  outweigh  the  nonlinearity. 

For  the  spatial  frequencies  in  the  interference  pattern,  this  condition 
prevails  over  almost  all  of  the  125-300p  region  in  which  the  last 
two  Brillouin  e-foldings  (and  hence  significant  backscatter  profile 
modifications)  are  taking  place.  The  small  decrease  in  widths  of  the  pump 
envelope  at  zl  {dQ  (B  =  5.6)  =  82  u,  compared  to  dQ(B  =  0)  =  85  u  from  Fig. 
and  the  backscatter  at  z2  (  (B  -  5.6)  =  75  V,  compared  to  d ^ (B  =  0) 

=  80  from  Fig.  10a)  shows  that  whole  beam  self  focusing  is  relatively 
unimportant  in  this  example.  However,  at  the  low  spatial  frequencies 
corresponding  to  the  such  widths,  expression  (31)  gives 0  sp'*'  K  ~  Wdy 
One  should  therefore  expect  whole  beam  effects  to  appear  if  either 
dQ(B  =  0)  or  dj(B  =  0)  become  significantly  smaller.  This  point  is 
illustrated  in  the  next  set  of  examples. 

Conditions  in  Fig.  11a  are  similar  to  those  in  Fig.  10a,  excepting 
that  the  pump  beams  are  now  focused  at  z2  rather  than  zl.  With  B  =  0 
and  fa  G  =  12,  this  configuration  retraces  as  well  as  the  earlier 
examples,  and  gives  a  correlation  H  =  .52.  The  total  width  of  the  back- 


10a}, 
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scatter  envelope  at  z2  is  now  only  60  p,  however,  and  a  deterioration 
of  the  retrace  due  to  whole  beam  self  focusing  soon  becomes  apparent  as  B 
Increases.  For  8  ■  2.77  and  in  G  ■  12  (Q»  2),  the  middle  peak  decreases 
by  25%,  and  the  one  on  the  right  by  50%.  For  B  =  5.17  and  £,n  G  =>  9*^ 

(Q  ■*  1),  the  retracing  has  almost  completely  disappeared,  as  seen 
in  Fig.  11b,  and  the  correlation  has  dropped  to  H  ®  .107.  The  whole 
beam  self  focusing  responsible  for  this  is  clearly  evident  at  the  z2 
plane. 

Finally,  Fig.  12  illustrates  the  importance  of  the  holographic 
mechanism  in  small  scale  ray  retracing.  Here,  the  inhomogeneous  plasma 
lies  between  zl  =  900p  and  z2  =  1200u,  so  that  the  pump  beams  remain 
almost  completely  separated  from  one  another.  Each  one  therefore 
generates  its  own  whole  beam  backscatter  from  the  random  noise  field 
in  its  vicinity,  as  seen  at  the  z2  plane.  Since  these  backscatter 
amplitudes  are  more  or  less  random,  and  their  phases  remain  uncorrelated, 
they  will  combine  to  produce  only  random  patterns  back  at  the  lens. 
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V.  SUMMARY  AND  CONCLUSIONS 


In  thi s  report ,  we  have  performed  numerical  studies  of  optical 
ray  retracing  in  stimulated  backscatter  from  laser-produced  plasmas. 
Ray  retracing  was  found  to  be  influenced  by  several  factors,  such  as 
the  optical  diffraction  lengths,  gain  inhomogeneity,  gain  narrowing, 
self  focusing,  and  the  size  and  location  of  the  focal  spot. 

For  whole  beam  retracing,  which  arises  from  the  high  axial  gain 
within  a  long  narrow  focal  waist,  we  found  the  following  results: 

(i)  In  simulations  where  a  half-blocked  beam  is  focused  near  the 
critical  surface  of  an  inhomogeneous  plasma,  only  a  small  portion 
(  5%)  of  the  backscatter  returns  to  the  lens  in  the  shadow  region. 
This  demonstrates  that  the  reflection  is  non-specular ,  and  is  in  good 
qualitative  agreement  with  a  recent  NRL  experiment.^  (ii)  Because 
of  gain  narrowing  effects  within  the  focal  waist,  the  spatial  profile 
of  the  backscatter  at  the  lens  tends  to  have  broader  wings  than  that 
of  the  pump  beam,  (iii)  The  backscatter  profile  at  the  lens  is 
relatively  insensitive  to  self  focusing,  provided  that  the  Brillouin 
gain  is  sufficiently  high  (e.g.  G>100).  At  lower  gain,  self 
focusing  can  cause  the  retracing  to  exhibit  a  pronounced  threshold 
with  respect  to  variations  of  density  or  pump  irradiance,  especially 
when  the  plasma  is  placed  just  beyond  the  focal  plane.  (iv)  When 
the  target  is  moved  far  out  of  focus,  so  that  the  beam  width  becomes 
comparable  to  the  length  of  the  gain  region,  the  retrace  effect 
disappears,  and  the  backscatter  profile  becomes  essentially  that  of 
ampl if  led  noise. 
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Small  scale  retracing  or  wavefront  reversal  occurs  when  the 

transverse  structure  of  the  pump  beam  creates  a  broad  interference 

pattern  in  the  plasma,  which  then  behaves  as  an  active  volume  hologram. 

According  to  this  model,  effective  wavefront  reversal  requires  that  the 

2 

gain  occur  slowly  over  typical  diffraction  lengths  1/k^^  • 

Excellent  agreement  was  found  between  a  simple  coupled  wave  analysis 
1 1  24 

of  the  model,  ’  and  the  numerical  solutions  of  an  idealized  test 
problem  where  two  truncated  plane  waves  formed  an  interference  pattern 
in  a  uniform  plasma. 

In  more  realistic  simulations,  we  chose  a  pump  profile  consisting 

of  three  isolated  hot  spots  at  the  lens,  producing  an  interference 

pattern  ofaslOOp  width  near  focus.  These  simulations  produced  the 

following  results:  (i)  wavefront  reversal  is  enhanced  if  the  plasma 

has  a  monotonical ly  decreasing  gain  coefficient  away  from  the  target. 

For  example,  a  plasma  density  variation  exp(-z/i.),  where  i:>>1/k0ej)  » 

produces  significantly  better  wavefront  reversal  than  a  uniform  slab 

of  thickness  Z  and  comparable  net  gain.  This  enhancement  arises  from 

the  better  spatial  filtering  action  that  occurs  in  the  lower  gain 

region  z  >Z.  (ii)  Gain  narrowing  effects,  which  were  present  in  all 

cases,  restrict  the  number  of  periods  within  the  backscatter  interference 
1 1  29 

pattern.  ’  At  the  lens  plane,  this  tends  to  eliminate  higher  spatial 

frequencies  and  broaden  the  hot  spots  in  the  backscatter.  (iii)  In  the 

case  of  a  plasma  dens  i  ty  ~  exp  (-z/Z),  self  focusing  appears  to  have  a 

2 

minimal  effect  on  the  wavefront  reversal,  provided  that  Z»  1/kQ6  D 
and  the  width  of  the  backscatter  Interference  pattern  is  large  in  com¬ 
parison  to  its  period.  If  this  width  decreases  below  a  critical  value, 


then  whole  beam  self  focusing  takes  over  and  destroys  the  wavefront 
reversal.  This  is  especially  likely  in  cases  where  gain  narrowing  effects 
are  prominent,  or  where  the  pump  beams  are  focused  a  few  hundred  microns 
in  front  of  the  target,  (iv)  If  the  target  is  moved  into  the  quasi  near 
field  of  the  lens,  so  that  the  hot  spots  no  longer  interfere,  the  wavefront 
reversal  disappears,  and  the  backscatter  profile  assumes  a  random  noise 
character.  The  formation  of  a  volume  hologram  is  therefore  essential  for 
wavefront  reversal. 

So  far,  this  work  has  been  restricted  to  the  strong  damping 
limit,  and  pump  depletion  has  been  ignored.  In  a  future  study,  we 
will  attempt  to  incorporate  phase  mismatch,  pump  depletion,  self  consistent 
ion  heating  and  ion  wave  saturation  into  the  code,  and  to  simulate 
actual  phase  conjugation  experiments.  As  suggested  in  Refs.  11  and  29, 
effects  such  as  pump  depletion  and  ion  wave  heating  are  likely  to 
improve  the  ray  retracing  by  reducing  spatial  gain  narrowing. 
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OPTIC  RAY  RETRACING 


Optical  ray  retracing  test  (Ref.  10).  No  backref lection  exposure  is  found  on  the 
blocked  half  of  the  beam,  indicating  that  ray  retracing  occurs. 


Fig.  4  -  Transverse  intensity  profiles  for  two  plane  hyper-Gaussian 
pump  waves  (6  =  o,  6  1  s  6  i.16  radian)  and  a  backscatter  "noise" 

wave  (6  =  0)° interacting0 within  a  uniform  100m  plasma  slab  center 

ed  at  z1  =  0.  In  this  and  all  subsequent  profiles,  the  solid  and 
dotted  lines  are  the  normalized  pump  and  backscatter  intensities 
respectively. 


POWER  GAIN  xio-4 


ANGULAR  DETUNING  01  -90  (RADIANS) 

Fig.  5  -  Tuning  curve  for  the  100V  slab,  showing  Brillouin 
gain  G  vs.  angular  detuning  9^  -  0  .  The  solid  line  is  the 
analytic  result  from  Refs.  11  or  24u (exact  for  infinite  plane 
waves),  and  the  dots  are  the  numerical  results  from  BOUNCE. 
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Fig.  6a-  Transverse  profiles  of  the  pump  radiation  and  backacatter 
from  a  300 P  plasma  of  density  .2  nc  exp  (-z/100)»  showing  partial 
retracing  of  a  half  blocked  pump  beam  at  the  lens.  The  stochastic 
noise  appears  as  random  dots  at  the  focal  plane  zl  ■  0.  CNegligible 
self  focusing  (B  =  0) ,  with  In  G  ■  113 
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Fig.  6b  -  Transverse  profiles  of  the  pump  radiation  and  backscatter 
from  a  300 y  plasma  of  density  .2  nc  exp  (-z/100),  showing  partial 
retracing  of  a  half  blocked  pump  beam  at  the  lens.  The  stochastic 
noise  appears  as  random  dots  at  the  focal  plane  zl  -  0.  [Strong  self 
focusing  (B  ■  7.6)  with  comparable  gain  £n  G  ■  8.93  The  density 
variations  n/n  are  given  by  Eq.  (15). 
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Fig.  7a  -  Transverse  profiles  for  a  uniform  400  plasma  located 
just  beyond  the  focal  plane  z2  ■  0.  CBelow  threshold  for  retrace 
(B  *  1.2,  in  G  =  1.3)  with  nQ  =  .01  ncJ 


39 


-.5  0  +.5 
x/f 

Fig.  7b  -  Transverse  profiles  for  a  uniform  400  plasma  located 
just  beyond  the  focal  plane  z2  =  0.  C Above  threshold  (B  ■  2.9, 
In  G  =  7.3)  with  the  same  incident  pump  Intensity,  but  Uq  ■  .0 
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Fig.  8  -  Transverse  profiles  for  a  uniform  100y  plasma  in  the  quasi 
near  field  of  the  lens.  The  whole  beam  ray  retracing  is  ineffective 
when  the  width  of  the  gain  region  becomes  comparable  to  its  length. 


Fig.  10a  -  Transverse  profiles  for  a  300w  plasma  of  density  .2  n^. 
exp  (-z/100),  showing  nearly  complete  retrace  of  the  three  pump 
beams  at  the  lens.  CNegllglble  self  focusing  (B  =  0)  with  In  G  * 
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Fig.  10b  -  Transverse  profiles  for  a  300  u  plasma  of  density  .2  n 
exp  (-z/100),  showing  nearly  complete  retrace  of  the  three  pump  beams 
at  the  lens.  C Strong  self  focusing  (B  ■  5.6)  with  In  G  «  103 
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Fig.  11a  *  Transverse  profiles  for  a  300  p  plasma  of  density 
.2  ng  exp  {-(z  +  300)/100}.  [Negligible  self  focusing  (B  *  0) 
with  J.nG  -  12] 


Fig.  lib  -  Transverse  profiles  for  a  300p  plasma  of  density 
.2  tic  exp  {-(a  +  300)/100),  showing  the  degradation  of  small 
acale  retrace  by  whole  beam  self  focusing.  (B  *  5.17  with 
InG  -  9.4) 
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Fig.  12  -  Transverse  profiles  for  a  300y  plas 
density  .2  n£  exp  {-(z  -  900)/100},  illustrat 
that  small  scale  retrace  disappears  when  the 
beams  do  not  overlap. 
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